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OO . 1 Introduction and summary 

cn 

O ■ This research is being done together with Georges Ripka and Bojan GolU. We will show 

. how to gauge effective quark models with non-local quark interactions. First, however, 

^ ! let me state the reasons why we want to use such models: 



- Non-local regulators arise naturally in several approaches to low-energy quark dy- 



D . namics, such as the instanton-liquid model or Schwinger-Dyson calculations 0, 

^ ' presented at this Workshop by Dubravko Klabucar. For the derivations and appli- 

cations of non-local quark models see, e.g., P,^,|5|J^J^,3,P,CT ■ Hence, we have to cope 
X ■ with non- localities from the outset. 

■ - Non-local interactions regularize the theory in such a way that the anomalies are 



preserved fgJTJ] and charges are properly quantized. Recall that with other methods. 



such as the proper-time regularization or the quark-loop momentum cut-off p!^JT3[ 
the preservation of the anomalies can only be achieved if the (finite) anomalous part 
of the action is left unregularized, and only the non-anomalous part is regularized. 
If both parts are regularized, anomalies are violated badly Jl^. We consider such 



division rather artificial and find it quite appealing that with non-local regulators 
both parts of the action are treated on equal footing. 

With non-local interactions the effective action is finite to all orders in the loop 
expansion. In particular, meson loops are finite and there is no need to introduce 



another cut-off, as was necessary in the case of local interactions |jT5|,|T6|,|T3 • As the 
result, the model has more predictive power. 
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As Bojan Golli, Georges Ripka and WB have shown |T8[, stable sohtons exist in 
a chiral quark model with non-local interactions without the extra constraint that 
forces the fields to lie on the chiral circle. Such a constraint is external to the known 
derivations of effective quark models and it is nice that we do not need it any more. 



What is the price to pay for all these nice features? 



- The calculations are more complicated — an extra integration over an energy variable 
has to be performed numerically. 

- Noether currents acquire non-local contributions. These are very much wanted since 
they make the Noether currents and anomalies conserved. However, there is an am- 
biguity involved. The transverse parts of currents are not fixed and their choice is 
part of the model building. Recall that this problem has been known for a long time 
in nuclear physics, where the transverse part of the meson-exchange currents is am- 
biguous. It is not possible to get rid of this ambiguity when one gauges non-local 
models. An ideal solution would be to first gauge the underlying theory {e.g. QCD 
with instantons) and then derive an effective gauged quark model. This has not been 
attempted so far and we need to deal with transverse currents which are not unique. 

- Non-local interactions modify the current-quark interaction vertex. In addition, con- 
tact terms (sea-gulls) are present in processes with more than one current. 

In this work we adopt the so-called P-exponent prescription ||19|,|9|jl0| for constructing 
Noether currents, with a focus on application to solitons. In particular, we show that: 



Noether currents corresponding to symmetries are conserved. In particular, the CVC 
and PCAC relations hold. 

Charges of the soliton (baryon number, isospin, gj^) do not depend on the choice of 
the path in the P-exponent, hence are unambiguously determined. They pick up a 
non-local piece, which is crucial for the charge quantization. 

Any n-point Green's function with external current momenta set to is independent 
of the path. An important example is the moment of inertia of hedgehog solitons, or 
transverse parts of vector correlators at zero momentum. 

Soliton radii, magnetic moments, in general form factors do depend on the path, 
hence are not uniquely determined. We show that this dependence is weak in the 
weak-nonlocality limit. 

A popular choice of the path in the P-exponent is the straight line P,|lO|- We give 
explicit expressions for evaluating Noether currents with this prescription which can 
be used in soliton calculations. 



2 The model 



We are concerned with the chiral quark model with non-local interactions, such as 
discussed in the talks by Georges Ripka and Bojan Golli. The Lagrangian is given by 
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L = y {dp)'ip{—p){p • J — m)ip{p) — J {dp){dp^)ip{—p)r{p)^^{p — p^)Tar{p')'ip{p') 

2 

j{dp)^\-p)^a{p). (1) 

For the purpose of this work we find more convenient to work in the momentum repre- 
sentation. We use the notation [dp] = J etc. The fields ^p describe the quark, m is 

the current quark mass, r{p) = r{p'^) is the regulator (local in the momentum space), 
and = / d^xe~'^^'^^"'{x) is the Fourier transform of the soliton field, which is 

local in the coordinate space. The index a = corresponds to the a field, with Fq = 1, 
and a = 1, 2, 3 denotes the pion, with F^ = i'y^Ta- We will also use the abbreviation 
$ = $'^Fa. The Euler-Lagrange equations have the form 



{p--f -m) %l){p) = J {dp')r{p)^''{p - p')Var{p')ip{p'), 
i>{-p) (p' ■ 7 - m) = J {dp)ilj{-p)r{p)^''{p - p)Tar{p), (2) 
'^"(^) = ~~^J (^^'OV'l-g - p'Viq + p')T^r{p')ij{p'). 

We note that by "unbosonizing" the model, i.e. by inserting the third of equations 
into Eq. (|I]), we recover the form with the quartic separable quark interaction given by 

L = j {dp)i}{-p){p ■ - m)i){p) + ^ J {dpi){dp2){dp^){dpi) x 

+P2 + + PA)i'{pi)r{pi)T ar{p2)'ip{p2)i'{p^)r{p^)T ar{p4)i){pA). (3) 

On the other hand, integrating out the quark fields from Eq. (|l]) leads to the bosonized 
(Euclidean) action, as used by George Ripka and Bojan GoUi in their talks: 

/($) = -Tr log idr-ia-V + (3m + /3r$r) + ArTr^^ (4) 

20"^ 

where Tr denote the full trace, i.e. functional as well as over color, fiavor, and Dirac 
space in the case of quarks. The forms (|l]), (|^) and (H) are fully equivalent. 

So much for the model. Its vacuum sector, spontaneous breaking of the chiral symmetry, 
the method of fixing the model parameters by fitting the pion mass and decay constant, 
and so on, are described, e.g., in Refs. |]9|JT0|JT8|1 . 



3 Gauge transformations 



The Noether construction of currents produces a contribution whenever a derivative 
acts on a field in the Lagrangian. In Eq. (|l]) the first term (the local term) involves 
one derivative, and results in the usual contribution to Noether currents. However, the 
interaction term with functions r{pi) may be viewed, in the coordinate representation, 
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as involving infinitely many derivatives acting on the quark field. This leads to com- 
plications. Below we show how to gauge the model in this case. Also, the presence of 
infinitely many derivatives does not allow for the canonical quantization of the quark 
field. We do not know how to quantize the quark field, and yet, as we will show further 
on, the charges, such as the baryon number, will be quantized. 

Let us consider the gauge transformations of the quark field: 

^(x) ^e-^^^-^-^^Vla;), (5) 

where for the cases of interest A° = \/Nc (baryon current), A'^ = r"/2 (isospin current), 
= 75''"'^/2 (axial current). The phases (t)a{x) parameterize the transformation. The 
local contribution to the Noether currents, i.e. the contribution coming from the first 
term in Eq. (|I]) is, of course, ja'^{x) = il){x)\°''')^il){x), or, in the momentum represen- 
tation, 

j,^'L(g) ^ j d'x e-^'^-^ji^'^^ix) = J {dp)^{-p)XaY^{p + q)- (6) 
With help of the equations of motion (0) we find that 

T /■ , + , \ rip' + q)\"-(3^{p' — p + q)r{p)— \ 

q-faiq)= idp){dp')^\-p'){ ^ (7) 

I r{p')f3^{p' — p + q)X''r{p — q) I 

Note that q ■ j]^{q) does not vanish in the non-local model! What is missing is the 
non-local part discussed below. 



4 P-exponents 



Now we are going to adopt a rather elegant way of gauging the non-local model. The 
P-exponent is defined as [p!9|,p|,p!0| 

W(x,i/) = Pe'/>''(^)''^^ (8) 

where A'^ = \°'A'^ is the gauge field (in general non-abelian), s parametrizes an (ar- 
bitrary) path from x to and P denotes the ordering along the path (needed only 
for non-abelian groups). Under the gauge transformation the field A'^ transforms as 
A'^ ^ A'^ + d'^(j)a- The following object transforms properly under the gauge transfor- 
mation: 

T = - J d^x J d^y J d^z tlj\x){x\r\z)W{x,z)(3^z)W{z,y){z\r\y)^{y). (9) 

Hence, when gauging the non-local model, we replace the interaction term in Eq. (|I]) 
with Eq. (I). 
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The non-local contribution to Noether currents, defined as ja'^^{(l) = ^r/(5y4^(— g)|^=o 
is equal to 

xti;\-h)rip) e-^'^-'ds^ A„/?$(/) + J e-^'^'^ds, /?$(/)A,| r(p')^(A;2). 

(10) 

Note that this expression is not unique, which is manifest in the freedom of choice of 
the path in the s integration. 



5 Two tricks 

In the following we shall use the following obvious formulas: 

h q-b 



q^" J e-^^ 'ds^ = J e-^^ 'diq ■ s) = te-"' '' - ie"'''", (11) 

a q-a 
b 

J ds^ = b^ - a^. (12) 

a 

The point here is that, clearly, the right-hand-sides of the above formulas carry no 



information on the choice of the path joining the end-points a and b. Formula (11 
appears when the longitudinal parts of currents are involved, and formula (|12|) occurs 
when the momentum of the current is zero. 



6 Conservation of currents 



With help of Eq. (^) we easily find the nonlocal contribution to the divergence of 
currents: 

[ r{p')l3^{p' - p + q)Xa[r{p - q) - r{p)] 

(13) 

Combining it with the local piece (^ yields the divergence of the total Noether current, 
q-Ja{q) = I {dp){dp')ip\-p')r{p')[\a, (3^p' - p + q)]r{p)^{p). (14) 
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For the baryon current this is immediately 0. For the isospin and axial currents we use 
the equations of motion for the $ fields to obtain 



d-Va{x) = -]^e'''''^,{x){i;r\x)V,{x\ri,) = -^e«''^$,(a;)$,(x) = 0, (15) 

d ■ Aa{x) = -^(^/'r|a;){r"75, ($o(a;) + m + i-fr,T^^b{x))}{x\r'ilj) = ma'^^a{x). (16) 

This verifies explicitly CVC and PCAC for the construction with P-exponents. Note 
that the the above expressions hold for any choice of the path. In other words, the 
longitudinal parts of vector and axial currents are fixed unambiguously. We add par- 
enthetically that this fact is related to the Ward-Takahashi identities, which hold in 
the non-local models. 



7 Charges 

Equation (^) simplifies greatly for the static case, g = 0. Then, through the use of Eq. 



(P^) we get 



j^'^^iO) = - d^x d^y d^z / {dh){dk2){dl){dp){dp') X (17) 



tlj^{-ki)r{p)Xa(3^{l)r{p')-iP{k2)-^e-''''-''+'''^-y+'^^^^^ 
< Pfj. 



^\-k,)r{p)l3<!>{l)Xar{p')^{k 



d 



-iki ■x+ik2-y+il-z+ip-{x—z)+ip'-{z—y) 



Next, we integrate by parts in the p and p' variables, denote r^{p) = dr^{p)/dp^^ and 
carry out the x, and z integrations. The result is 

r{p)l3^{p - p')\ar^,{p') 

When the quark fields are integrated out, the following formula for the full currents at 
g = holds: 

p[—io ■ 7 + r<I>r) 
For the charges, Qa = ja{0), we have 

= -Tr ^, \ ^ , iPloXa + roXaP^r + r/^^A^ro) . (20) 
p{—io ■ 7 + r<I>r) 

As advocated earlier, ja{0), or the charges do not depend on the path in the P- 
exponents. Certainly, this should make us happy. In local field theories the charge 
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is fixed by quantization. Here we can see a similar feature, although we have not quan- 
tized the quark field. 



8 Baryon number of the soliton 



We shall now examine the issue of the baryon number in some greater detail. For 
stationary solutions (such as solitons) the fields $ are time-independent. In this case 
(we pass to Euclidean space in this section) 



Qa = -E/ 



dv 



27r iu + Sk (v) 



{k, u\{Xa- iroXaP^r - ir/3$A„ro) \k,u), 



(21) 



where we have used the spectral representation of the energy-dependent Dirac Hamil- 
tonian 



— * 

h{i') = —ia ■V + m + /3r(i^)$r(i/), 
with the energy-dependent spectrum: 

h{i/)\k,i') = Ekiv) \k,i/). 
In particular, for the baryon number we obtain 

dv 1 



(22) 



(23) 



B 



1 



-oo 
oo 



27r iu + £k (v) 

dv 1 — i{k, u\hQ\k, u) 
27T if + Sk {y) 



[1 - i{k, i/|/3(ro$r r$ro) v)] 



" A, 1 - 

27Ti V — isk (v) 



(24) 



where ho = dh{u)/dv = /3(ro$r -|- r$ro), and the Feynman-Hellman theorem has been 
used in the last equality. 



Suppose we have a pole in the quark propagator at v = Vq, i.e. 



(25) 



{k labels quantum numbers relevant for the soliton, such as grand spin, parity, and 
radial number). Expanding the denominator around v = Vq we find 



1^0 + {v - Vo) - iSk iyo) - i 
dsk (v 



dsk (f ) 



[v - Vo) [I - i 



dv 



dv 
+ ... 



{V - Vq) + 



(26) 



v=vo 



V=VO'' 
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Fig. 1. The contour in the integral over v encirchng a positive-energy valence state. 

We notice that, quite remarkably, the numerator in (|2^) is such, that the residue of 
any pole is equal to unity. This means, that the baryon number is properly quantized 
in the model. Therefore, we have achieved the baryon quantization without quantizing 
the quark field! This is a very important feature, which brings the model close to the 
particle-hole interpretation: by occupying (iVc times) a pole of the quark propagator we 
raise the baryon number by one unit. We achieve this by distorting the contour in the 
V integration such that it encircles the occupied valence states (this is only necessary 
if these states lie above the real axis). See Fig. 1. 

It can also be shown straightforwardly that the energy of the stationary system equals 
to 

, J 2m V — le-k (v) 

hence occupying {Nc times) a state corresponding to a pole at v = vo brings the energy 
NcEk (vq) . The Dirac sea contribution is obtained from expression analogous to (p4D, 
but with the contour undistorted. 



Above we have said "c/ose to the particle-hole interpretation" for the following reason. 
Unlike the usual many-body problem, where all poles of the quark propagator lie on the 
imaginary axis (recall we live in the Euclidean space), in the non-local case there are in 
general many poles in the complex plane. In fact, this feature is also present in certain 
local quark models |]2D|. These poles are induced by the presence of the regulator. 



Since they are complex, they do form asymptotic states. In fact, it is possible to choose 
the regulator in such a way, that in the vacuum there are no poles on the imaginary 
axis (no "physical" poles), which is sometimes referred to as "analytic confinement". 
Remarkably, the hedgehog soliton fields generate a valence state on the imaginary axis 
with a low eigenvalue (Bojan GoUi's talk). It is therefore natural to occupy this state, 
thus providing the soliton a unit of baryon number. 

The contribution of the Dirac sea levels to observables {e.g. as in (|27|)) is obtained by 
carrying numerically the integration over the v variable along the real axis. 
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We also remark that expressions analogous to (^4]) hold for any "good" quantum num- 
ber, i.e. when [/i (i/) , Aq] = 0. For instance, in the Friedberg-Lee-like model, where 
$ = 0" (no pion field present), the isospin is a good quantum number, and we have, in 
analogy to (|2^) , 



1^ 7|^; (l-^^)(Mr3|fc,.) ^ (Freedberg -Lee) (28) 



2^ j 2711 V - lEk [V) 

— oo 

In hedgehog models isospin is not a good quantum number and it has to be restored 



by a suitable projection method, e.g. by cranking ||21| , |22 



9 gA 



Another important quantity is the axial charge of the nucleon, gA, evaluated in hedge- 
hog models from the expectation value of the z component of the axial current. In fact, 
we see that also the space components of current do not depend on the path at g = 0. 
Hence gA is path- independent. An explicit expression can be immediately derived from 
(pj|), using the method of Ref. . 



10 Moment of inertia 



It is well-known that hedgehog solitons break the spin and isospin, which are restored 
by a suitable projection method. In the cranking method 0j2^ the basic dynamical 
quantity is the moment of inertia ^. It is obtained by adiabatically rotating the soliton, 

^(x) ^ e-*"''^"V(a;), (29) 

where Qa is the (adiabatically small) velocity of rotation in the isospin space, and t is 
the time. The moment if inertia is obtained by performing the transformation (|2^) in 
the action and then identifying the coefficient of fl^. We notice that Eq. (|29D is a special 
case of gauge transformation (^, with the vector potential equal to = {fla, 0, 0, 0). 
Therefore we should apply the prescription (|^), and then evaluate 6 from the formula 
6"^! /6flaS^b = ^""^0. Applying the same techniques as used in the previous sections 
leads to the following expression for the moment of inertia in models with non-local 
regulators (Euclidean notation): 



X 



Q__\y- f {k, y\ (rg - iT^(3ro^r - z/^r^rpTg) |/, u) 
2f^J2TT iv^ek{y) 
(/, v\ (rg - irg/jro^r - i/jr^rpTs) \k, v) 
iv + Eiiv) 

^ r dv {k, u\(3 Uroo^r + rgro^roTg + ir^roo) \k, v) 
-Y r-. , 30 
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where roo = d?r/dv'^. The first term in Eq. (|30D is dispersive, it is involves to quark 
propagators. The second one is the contact term (sea-gull), with one quark propagator 
looped around. In the local limit Eq. (|30| ) reduces formally to the usual Inglis formula 

^ p,h 

In the presence of valence states we proceed as in the case of the baryon number or the 
soliton energy, and decompose the total moment of inertia into the valence and Dirac 
sea parts. The valence part is obtained by explicitly occupying the state satisfying Eq. 



Notice that 6 in Eq. (pO]) is path-independent. In fact, it is clear from the derivation 
that any n-point Green's function with vanishing momenta on the external current 
lines does not depend on the choice of the path. This is because the differentiation of 
the action with respect to the A field at g = brings down, according to Eq. (12), 



the factor of {y — x)^, where y and x are end-points of the line. Then A is set to zero. 
Obviously, no information of the choice of the path is left. 



11 Straight-line paths 



The quantities discussed above did not depend on the choice of the path. This is not true 
for other physical quantities, such as form factors and magnetic moments of baryons, 
or transverse vector correlators. The simplest choice of the path in the P-exponent is 
just the straight line, as used in Refs. HJl^l- One parameterizes = + a{z^ — x^). 
Then 



and the repetition of the steps of Eq. (|T^-[TB|) leads to 



(32) 



(dp) (dp') / daip\—p) X 



r^{p + aq)\aP^{p - p' + q)r{p') + 



r{p)/3^{p -p' + q)Xarf^ip' - aq) 
In the coordinate representation we can write equivalently 



(33) 



ja^^{x) = — J da J dz dx'6{x — x' — a{z — x')) x 



\z){z\r^\x'){x'\Xa/3^r\ip) + {'ip\r/3^Xa\x'){x'\rf,\z){z 



(34) 



10 



Expression (| 



,0) can be used to calculate the non-local contribution to currents. 



12 Form factors 



Rewriting the general expression ( p^ ) for \a = l/N^ we find the following expression 
for the non-local contribution to the Fourier-transformed baryon density: 



1 



ro{p + aq)^{p — p' + q)r{p') + 
r{p)^{p — p' + q)ro{p' — aq) 



We can now pass to the Breit frame (go = 0) and expand Eq. (|35| ) at g = 0. The term 
quadratic in q is related to the non-local contribution to the baryon mean squared 
radius, which equals to 



2\NL 



1 ^ f 

/ {dp){dp')H-p) X 



N, 



c i=l 



[ro{p)^ii{p - p') + roi{p)^i{p - p') + \rQii{p)^{j) - p')'^r{p') + 
r{p) [^u{p - p')ro{p') - - p')roi{p') + |<l>(p - p')roii{p')} 



(36) 



The expressions for the isoscalar and isovector magnetic moment involve cranking and 
are somewhat more complicated. They can be obtained using the methods described 
above along the lines of Ref. ||22|| . 



13 The weak-non-Iocality limit 



It is instructive to have a closer look at Eq. (|36|). Suppose the soliton has a typical 
size R ~ l/As', and the regulator has a momentum scale A above which the momenta 
are suppressed, e.g. r{p'^) = exp(— p^/A^). Derivatives with respect to momenta, such 
as appear in Eqs. (^), bring down a factor of the inverse scale squared, e.g. r^{p'^) = 
— exp(— p^/A^)/A^, and similarly for the soliton profile, where ^j{P — P') = 2{P — 
P')j/A| $'(P — P'). In the weak-non-locality limit, i.e. when A S> As, the terms with 
derivatives of the regulators are suppressed relative to the terms with derivatives of $. 
For instance, in Eq. ( ^6]) we must keep only the terms with and can neglect the 
remaining pieces. In the coordinate representation this is equivalent to using, instead 
of (^41), the following expression for the non-local current: 

^■M,NL(^) = - mx){x\r^Xaf3<!>r\^) + (V^|r/5$A,r^|x)(x|^)] . (37) 

One can formally pass from Eq. (0) to Eq. ( p7D by commuting the and |a;)(x| 
operators, which is allowed in the weak-non-locality limit Q. 
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For the solitons shown by Bojan Gohi A ~ IGeV and A5 ~ .25GeV, hence A|/A^ ~ 
0.06, hence we seem to be very close to the week-non-locahty hmit. This is fortunate, 
since then the observables such as the baryon radius, etc., do not depend strongly on 
the choice of the path in the P-exponent. 

Our research on solitons with non-local regulators is under way and we hope to be able 
present further exciting results shortly. 
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